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THEORETICAL ERROR ANALYSIS OF A DOPPLER RANGE-RATE 


AND PHASE -MODULATED RANGE TRACKING SYSTEM 
By Kenneth J. Bures and Gerald L. Smith 
Ames Research Center 


SUMMARY 


A detailed error analysis is presented for a particular type of range and 
range -rate tracking system. Results are given in the form of analytical 
expressions for the variances of the various types of errors that affect 
tracking data. The expressions are written in terms of system parameters so 
that they can he applied to a wide range of specific systems and tracking sit- 
uations. The majority of the errors considered are random, but some are bias- 
like and others are deterministic. The equations given constitute an error 
model from which data weighting factors can be obtained for use in data 
processing for the construction of the "best estimate" of the trajectory of a 
space vehicle. 


INTRODUCTION 


When data from a tracking system are used to estimate the trajectory of a 
space vehicle, whether by means of least squares, Kalman filter, or some other 
method, it is necessary to apply appropriate weighting to each data point to 
obtain a "best estimate." Since these weights should be proportional to the 
accuracy of the data, one must have a good idea of the errors that contribute 
to data inaccuracies. The method of obtaining this information is to perform 
a detailed error analysis of the equipment, that is, to identify as many of 
the error sources as possible and, from the known configuration of the system, 
determine the effects errors from such sources should have upon the data 
generated in the tracking operation. 

This report is concerned with developing, in the above manner, a detailed 
error model for a hypothetical interrogator-transponder-type sinusoidally 
phase -modulated range and Doppler range-rate tracking system. The system is 
assumed to consist of an interrogator aboard a. satellite orbiting the moon and 
a beacon transponder located on the lunar surface. However, the results are 
equally applicable to the more conventional arrangement of an earth-based 
interrogator and a space vehicle transponder. 

Sinusoidally phase -modulated range and Doppler range-rate systems of the 
type described here presently find applications in a space vehicle tracking 
system (refs. 1 and 2), in the SECOR and ANNA geodetic satellites, in the 
SHIRAN airborne electronic surveying system (ref. 3 ), and in a zero-zero air- 
craft landing system (ref. 4 ) . Despite their wide use, it appears that no 
complete, detailed error analysis has been performed on such range and 



range-rate tracking systems. A thorough search of the literature revealed 
only two previous analyses (refs. 5 and 6 ), both concerned with the accuracy 
of the range -rate measurement only. This report considers those errors in 
both the range and range rate which are listed below: 

Range -Rate Errors 

1 . Master oscillator short-term frequency instability 

2 . Master oscillator long-term frequency instability 

3 . Quantization error 

k. Phase -locked loop (PLL) error 

5. Uncertainty in Doppler averaging time 

6 . Uncertainty in the speed of light 

7. Averaging time error 

8 . Finite propagation time error 

Range Errors 

l . Master oscillator short-term frequency instability 

2. Master oscillator long-term frequency instability 

3 . Quantization error 

4. Phase -locked loop error 

5. Phase detector error 

6 . Calibration drift errors 

7. Uncertainty in the speed of light 

8 . Finite propagation time error 

9 . Doppler shift in the ranging tones 

Some of these errors are random ( i.e. , uncorrelated from one data point 
to the next), some are biases or bias-like, and some are deterministic. In 
the report these characteristics are pointed out, and formulas for the rms 
values of each error are developed. 
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PRINCIPAL SYMBOLS 


B 

c 

f(t) 

f c 

f L 

f Ld 

f t- 

f td 

G r 

Gb 

Gcp(co) 

k 

L(w) 

m 

N 

% 

N y 

P 

P r 

Pt 

Q 

r 

r B 

R 


noise equivalent bandwidth of master oscillator phase noise 
velocity of light 

random frequency component of master oscillator 

range counter pulse generator frequency 

transponder offset frequency 

one-way Doppler shift of f^ 

bias frequency added to 2f^ d 

interrogator transmitter carrier frequency 

one-way Doppler shift of f^ 

receiver antenna power gain 

transmitter antenna power gain 

power spectral density of phase fluctuation 

Boltzmann's constant 

second -order FLL transfer function 

frequency multiplication factor 

final count of pulse counter 

power spectral density to input of second FLL 

power spectral density to input of first FLL 

oscillator power 

received power (at antenna) 

radiated (transmitted) power 

oscillator selectivity 

radius vector to satellite 

radius vector to beacon 

true range 



R m 

ft 

R 

S L 

Ss 

t 


ti 


t 2 

T 

T c 


• • 

X,Y 


Y a ^Y a 

M # * 

X,Y 

Aa 

P 

e 

i 

e 

x 

v 


P 



measured range 

estimate of the true range 

autocorrelation function 

long-term rms fractional frequency stability, Af/f 
short-term rms fractional frequency stability 
time 

beginning of Doppler counting interval 
end of Doppler counting interval 
Doppler averaging time 
oscillator coherence time, 2/$ 
temperature, °K 
vector components of R 
vector components of R a 
vector components of fi 

elapsed phase angle of carrier during time T 
"beacon latitude 

phase error in quantization process 

PLL damping ratio 

true anomaly of satellite orbit 

beacon longitude 

wavelength of ranging tone 

2B 

component of R in orbit plane 
standard deviation of calibration drift error 
standard deviation of phase fluctuation 
standard deviation of phase detector phase error 
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°R 

°R 

T 

Cp(t) 

q5 d (t) 

^md 

$ 

$ 

*n 


w d 

co n 


standard deviation of range 

standard deviation of range rate 

two-way propagation time, approximated by 2R/c 

phase angle of noise component of oscillator output 

phase angle of Doppler shift frequency 

phase angle between received and reference ranging tone 

phase angle of a ranging tone with a Doppler shift 

white power spectral density of f(t) 

normalized power spectral density of additive white noise at the 
input of the first PLL 

normalized power spectral density of additive white noise at the 
input of the second PLL 

PLL damped natural frequency, - \ z 

PLL undamped natural frequency 



RANGE RATE 


Description of the Range -Rate System 

Figure 1 is a block diagram of the range -rate measuring system and 
figure 2 is a block diagram of the transponder. 

The satellite interrogator transmits a signal with carrier frequency ft* 
The signal received by the beacon transponder has a frequency of ft + ft d , 
where ftd is 'the one-way Doppler shift due to the relative motion of the 
beacon and the satellite. The received signal is mixed (fig. 2) with a local 


oscillator signal of frequency (m-l)fp 



Figure 1. - Range rate system. 



Figure 2.- Transponder. 


, where fp is the local oscillator 
frequency and m is a frequency multi- 
plication factor, to give an IF sig- 
nal with frequency f t + f t(i - (m-l)fp. 
The IF signal is amplified and then 
mixed with another local oscillator 
signal of frequency mfp to give a 
signal with frequency ft + ft d + fp. 
This signal (an amplified copy of the 
input signal shifted in frequency by 
fp) is further amplified and retrans - , 

mitted. In addition, a carrier with 
frequency mfp is transmitted. 

Because of the relative velocity 
of the satellite and beacon, the 
transponder signals are received at 
the satellite with frequencies 
ft + fL + 2f td + fpd and m(fL + fLd) . 
respectively, where fp<i is the one- 
way Doppler shift in the transponder 
local oscillator signal. The received 
signal is mixed with the interrogator 
transmitter signal producing sum and 
difference frequencies fj + 2f+^ + ft .a, 
2f t + f L + 2ftd + fLd. and 
m(fp J + fp^) ± ft (see fig. l) . Number 1 
IF amplifier passes only the signal 
with frequency fp + fp^ + 2fp d , and 
number 2 IF amplifier passes only the 
m(fp + fp d ) - ft signal. Through a 
series of signal manipulations, repre- 
sented in figure 1 as simply PLL 
(phase -locked loop) but described in 
detail in reference 1, signals with 
frequencies fp + fpd + 2ftd and 
f L + fLd " fo. where f Q is a known 
bias frequency, are obtained. These 
are mixed together and low-pass 
filtered to yield a signal with 
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frequency 2fta + f Q , which is just the two-way Doppler shift of the 
interrogator carrier plus a known bias frequency. 


This frequency is measured by counting the number of 'zero crossings of 
the Doppler waveform in a given interval of time determined by the master 
oscillator. The count is related to the two-way Doppler shift; therefore, the 
range rate can be obtained in the computer by the relationship between range 
rate and Doppler frequency shift. That relationship, derived in appendix A, 
is 


cf. 


R = 


td 


ft 


( 1 ) 


where R is the range rate and the two-way Doppler frequency shift is 2ftd* 


Stochastic Errors 

Master oscillator short-term 1 frequency instability . - If the transmitter 
signal of frequency ft drifts by an amount Aft during the roundtrip propa- 
gation time, then on reception the apparent Doppler shift will be 2ft<i + Aft. 
The Aft is interpreted as part of the two-way Doppler shift and hence an 
error in the range-rate measurement occurs. 

Suppose the transmitted signal has a frequency ft + f(t) where f(t) is 
a zero mean time -stationary random frequency fluctuation. Then the signal 
received by the interrogator is of frequency ft + fL + 2ftd + fLd + f(t-r), 
where t is the two-way propagation time. The output of the first mixer is 
equal to the difference frequency f^ + fpd + ^fta + f(t-r) - f(t), After 
several stages of signal processing as described previously, the offset fre- 
quency fj, and its associated Doppler shift fpa are removed and a bias fre- 
quency f Q is added. The bias frequency serves as a convenience in the cycle 
counting process and is removed in the computer. Therefore, the following 
analysis is not affected if this term is neglected and the input to the 
counter is considered to have a frequency + f(t-T) - f(t). Associated 

with this frequency . is a phase angle cp^ + cp(t-r) - cp(t). 

The cycle counting process, in effect, measures the elapsed phase angle 
Aa of the input signal during the time interval t 2 - ti = T, where ti and 
t 2 are the times at the beginning and end of the counting interval. The 
average frequency is 


| [<Pd(t 2 ) + <p(t 2 - t 2 ) - <p(t 2 ) - <Pd(ti) - cp(ti - T X ) + cp(ti)J (2) 


1 In this report, "short-term" means time durations comparable to the 
two-way propagation time t. 
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This average frequency is interpreted as the Doppler shift frequency and, 
therefore, can he substituted for in equation (l) to yield the apparent 

average range rate over period T: 


_ c Acp d 
" 2w t T 


2(*>tT 


[cp(t 2 ) - Cp(t 2 - T 2 ) - cp(tl) + Cp(ti - Ti)J 


(3) 


The second term of equation ( 3 ) is the range-rate error due to oscillator 
instability. The rms value of this error is 


= 


R 2w t T 


(0>(t 2 ) - 9(tx)J - cp(t 2 - t 2 ) - Cp(tl - Ti) Y) 


(M 


where the bar means "average." Equation (4) can be expanded within the 
brackets to yield 


ffA = 


2od-fcT 


-j 9 (t 2 ) - cp(ti) - 2^cp(t 2 ) - cp(ti)][cp(t 2 - t 2 ) -cp(ti- ti)J 


[cp(t 2 - T 2 ) - cp(tl - T l)J^ 


(5) 


The averaging time T is normally kept small enough so that ti = r 2 = r. 
Using this approximation and noting that t 2 = ti + T we see that 


and 


cp(t 2 ) - cp(ti) = cp(ti +-T) - cp(ti) (6) 

cp(t 2 - t 2 ) - cp(ti - ti) = cpjjtx - t) + TJ - cp(tx - t) (7) 


Equation (6) represents the phase difference in the Doppler signal between 
time tx and t 2 , an interval of T seconds. Equation (7) also represents the 
phase difference which occurs in an interval of T seconds; but in this case, 
the interval T is shifted by r seconds with respect to that of equa- 
tion (6). Therefore, the first and last terms in the brackets in equation (5) 
are both equal to ft( 0) , the autocorrelation function of the phase difference 
with a correlation time of zero seconds; and the middle term in the brackets 
in equation ( 5 ) is equal to 2ft(r) , twice the autocorrelation function of the 
phase difference with a correlation time of r seconds. Therefore, 
equation ( 5 ) can be written as 
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I 


a R = 


= 7 =£ — f"ft(0) - R(t)"| 1/2 

v2oJ t T L J 


( 8 ) 


The autocorrelation function ft( t) must now be calculated. 

Frequency is defined as the time derivative of phase; therefore, the 
phase difference Acp of equations (6) and (7) can be written in terms of 
frequency as 


Ap = 


f(t)dt 


(9) 


t-T 


Integrating by parts yields the Fourier transform of equation (9). From this 
the power spectral density of the phase difference is found in terms of the 
power spectral density <5(w) of the frequency fluctuation f(t). Finally, the 
inverse Fourier transform yields the autocorrelation function Sl(r). (For 
details, the reader is referred to ref. 5») When the frequency fluctuation is 


white noise with a one-sided power spectral density 
function of the phase difference is 

$ , the 

autocorrelation 

R,(t) =< 

r 

$T 

2 

(1 - 0 < H 

< T 

(10a) 


0 

T < | T | 

< 00 

(10b) 


The autocorrelation function ft(0) is obtained by letting r = O' in equa- 
tion (10). Substituting R,( t) and R,(0) into equation (8) and letting 
2/$ = T c , the oscillator coherence time, 2 results in an rms range-rate error 
equal to 


_c_/2R_y /2 
s/iwtT \ cT <v 

c f . t_y ' 1 2 

s/2u) t T vV 


0 < R < 



< 00 


(11a) 


(lib) 




The white noise assumption is, of course, a simplification. Baghdady, 
Lincoln, and Nelin (ref. 8) have shown theoretically that a more realistic 
power spectral density for f(t) is 

2 Coherence time (ref. 7) is defined as the time required for the rms 
value of the phase difference between a noisy oscillator and a perfect 
oscillator to reach one radian. 
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( 12 ) 


$(w) = ^ + b + c w 2 

except for High and very low frequencies. The first term of equation (12) is 
flicker noise , a low-frequency process caused by current and voltage fluctua- 
tions in electron tubes and transistors, by mechanical vibrations, and by tem- 
perature variations. The second term of equation (12) is white noise due to 
thermal noise in the oscillator. Finally, the third term of equation (12) is 
due to noise associated with certain accessory circuits of the oscillator. 
Usually, the response of the buffer amplifier following the oscillator falls 
off faster than l/w 2 for frequencies greater than half the amplifier band- 
width, so that #(w) gradually approaches zero as w approaches infinity. 
Equation (12) predicts an infinite power spectral density at w = 0; however, 
in practice, $(u>) levels out rapidly and may even decrease as w approaches 
zero (ref. 8). Therefore, $>(u)) is not as dependent on frequency as equa- 
tion (12) would indicate. This is substantiated experimentally; in fact, 
Vessot, Mueller, and Vanier (ref. 9) have found that thermal noise is the 
major factor limiting short-term oscillator stability, particularly with 
atomic standards . 

Equation (ll) can be arranged to show the dependence of °R on the 
oscillator short-term stability factor S s and on the range R by eliminat- 
ing the coherence time T c . Substituting Edson's expression of coherence time 
(ref. 7) in terms of oscillator parameters, namely. 


T c 


2PQ 2 

KT°u>t 2 


(13) 


into his expression for rms fractional phase deviation. 


Ap 


-I, 


KT 


2TPQ 2 


(14). 


results in an expression for the fractional phase (and, hence, frequency) 
deviation in terms of coherence time: 


Aj> _ Af _ „ _ 1 

cp " f s 


(15) 


Finally, substituting equation (15) into(ll) and letting r = 2R/c gives the 
desired result: 
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o < R < Y ( l6a) 


7T < R < 0° (16b) 


Master oscillator long-term frequency instability . - In addition to the 
white -noise -type, short-term frequency fluctuations discussed in the previous 
section, a bias -type error is introduced when the Doppler frequency shift is 
converted into range rate. Equation (l) shows that the oscillator frequency 
ft must be known exactly if the range rate is to be determined correctly. 
However, since ft is not precisely known, an error occurs. 

The range-rate error is found by differentiating equation (l) with 
respect to ft and substituting the derivative into the total differential: 


/Rc\ 
\T J 


1/2 




CSg 


AR-|BjAf t 


(17) 


The result is 


AR = 


cf td Af t 
It It 



(18) 


If Aft/ft is interpreted as the rms uncertainty in ft, that is, the rms 
long-term fractional frequency instability' S-^, then the rms range-rate 
error is 

a E = rs l (19) 


Quantization error .- The process 
Doppler frequency in a given time T 



Figure 3» _ Quantization error. 


of counting the number of half-cycles of 
can be viewed as an analog-to -digital 
conversion of a phase angle. .However, 
the conversion is not exact. The 
counter % as a resolution of only one 
half -cycle since it counts zero cross- 
ings; therefore, the half -cycle count 
N does not include the portion e 
shown in figure 3* Because the phase 
angle e is not included in the count 
N, the measured frequency is less than 
the true frequency. All values of e 
within the bounds 0 < e < jt are 
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possible and equally likely. The maximum error, e mfl Y = it, results in the 
max imum count error ANmay of -1. The count N is equal to 


N = 2T(f 0 + 2ftd) 


( 20 ) 


The maximum error in Doppler frequency Afta mpy can be obtained by evaluating 
the total differential 


2 Af tdmax 


2 * f td 
dN 


ANjjjax 


AN, 


max 


2T 


2T 


( 21 ) 


where df^ d /dN is obtained from equation (20). The maximum range-rate error 
AR mav can be obtained by differentiating equation (l) with respect to f td 
and substituting the derivative into the total differential: 

• Sr 

AEmax = df td ^td^ ( 22 ) 

Then, substituting equation (2l) into (22) yields an expression for the 
maximum range -rate error: 

A^max = 4f^r ( 2 3) 


Since all e, 0 < e < it, are equally likely, the probability density 
function of e is uniform. Therefore, AR also has a uniform probability 
density function with probability l/AR ^y over the interval 0 < AR < j AR rnfly . 
It can be shown that for this rectangular distribution the mean is A Rmay /2 
and the standard deviation is AR mfl y/2 \ [ 3 . 


The error just considered also appears at the beginning of the counting 
period. In this case, the error is always such that the frequency reading is 
greater than the true value and the maximum count error is now A Nma.y = + 1 . 
The analysis follows exactly as above, with the results that the mean is 
-ARjjjax / 2 and the standard deviation is ARmax / 2 */ 3» 


The quantization errors at the beginning and end of the counting period 
are independent; therefore, their variances add. Thus, the overall 
quantization error variance is 
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(24) 


Substituting equation ( 23 ) into (24) gives the rms error in range rate due 
to quantization: 


12 



c 


(25) 


04, = 


R 4\f6f t T 


Note that the rms range -rate error is not a function of range or range rate; 
it depends entirely on the digital -to -analog conversion process. The mean 
quantization error is zero. 

Phase-locked loop error .- This section considers the effect of the phase- 
locked loop (PLL) on the range -rate measurement. The phase -locked loop intro- 
duces error in the range-rate measurement in two ways: (l) by the inability 

of the loop to track the Doppler signal exactly and (2) by the additive noise 
which corrupts the input signal and is subsequently processed with the signal 
by the phase-locked loop. 

The first error can usually be made negligible by proper design of the 
loop filter. With a high-order loop filter, the error could theoretically be 
exactly compensated since the Doppler shift is a more or less deterministic 
function. 

The effect of additive receiver noise cannot be neglected. There are two 
types of additive noise interferences: thermal (white) noise in the receiver 

components and external noise sources. In conventional communications systems 
operating at frequencies above 300 MHz the most important source of noise is 
the input section of the receiver (electron tubes, transistors, and crystal 
mixers) . However, with receivers employing extremely low noise devices 
(cooled crystal mixers and masers), the noise attributed to external sources 
is the limiting factor in receiver performance. The sun and the Milky Way are 
the two principal external noise sources; 3 of these two, solar noise is sev- 
eral orders of magnitude greater than galactic noise. The solar noise con- 
sists essentially of a slowly varying component (sunspot cycle) superimposed 
on a stationary white noise component. The galactic noise spectrum is not 
white, but over the relatively narrow bandwidth of the receiver the spectrum 
is essentially constant (see ref. 10 for further information). 

Thus, in the following analysis the noise is assumed to be gaussian white 
noise, with a constant power spectral density Ny. The value of Ny depends 
to some extent on the carrier frequency. However, any accurate figures on the 
noise power must be obtained by actual measurement, or at least by analysis of 
the particular equipment and mission under consideration. The values given in 
this paper are typical. 

It can be shown by letting cp(ti - ti) and cp(t^ - t 2 ) in equation (4) 
equal zero that the rms error in range rate due to phase fluctuations 
(rather than fluctuations in phase difference as in eq. (4)) is 

°R = 2^T + T) 2 - 2cp(t)cp( t + T) + cp(t) 2 ] 1/2 (2 6 ) 


3 For an earth-based receiver the earth's atmosphere would also be a noise 
source. 


13 



where the bar means average, and cp(t) represents the phase noise due to 
additive white noise. In terms of autocorrelation functions, equation (26) 
becomes 


Pa = r=r~ U(0) - R(T )] 1/2 
K v 2u t T 


(27) 


The autocorrelation function R(T) is found by the well-known technique of 
taking the inverse Fourier transform of the power spectral density of the PLL 
output. This is written in terms of the closed-loop transfer function L(w) 
and the phase noise normalized power spectral density $ n as 



|L(w) | 2 $ n (u))e^ a) ^' dw 


(28) 


Input 

signal 



(a) Block diagram. 


A block diagram of the phase -locked 
loop is shown in figure 4(a) and a lin- 
earized model is shown in figure 4(b). 
For the linear proportional-plus - 
integrator filter F(s) shown in fig- 
ure 4(c), the closed -loop transfer 
function L(w) of the PLL is 



(h) Linear model. 


L(w) 


V 5 + j2Q n u 
u^ 2 + j2f;u) n co - w 2 


(29) 


where w n is the undamped natural 
frequency and £ , the damping ratio of 
the system. Hence , 


rA/W— I 


-AAAAr 


|L(«) 


i+hff 


F(s) - 




GO 4 - 


Amplifier 


1)1 rr- ) +i 


(30) 


o — — — ; — — — o 

(c) Active linear filter , F(s). 
Figure 4. - Phase -Locked Loop (PLL). 


substituting equation ( 30 ) into ( 28 ) 
and integrating using residues 
(ref. ll ) 4 gives 


Reference 5 gives ft( T) for a second -order PLL containing a passive 
linear filter. Here R(t) is presented for a second-order PLL with an 
active filter because the active filter provides better tracking performance 
(p. 29 , ref. 12 ). 









®n“n " W r£ T| 


R(T) = e 


(1 “ ^ 2 )sin w d T + 


1 + 


4( 2 ' 

cos oo^T 


(31) 


where w d is the damped natural frequency w n "Jl - £ 2 . The autocorrelation 
function ft(0) is found by letting T = 0 in equation (3l) 


R(0) 


^n^n 1 + 4-t 2 

8 i 


(32) 


Equations (3l) and (32) are substituted into equation (27) to give the rms 
range-rate error due to additive white receiver noise: 


R 


C fl + 4t, 2 

8jtf t T \ i 


-^Tr 


“n 

w d 


(l - 4£ 2 )sin w d T + 


1 + 4< ; 

i 


cos 



(33) 


The normalized power 
as 


n 


spectral density of the additive white noise is defined 


N 

Power spectral density of noise _ V 
Average received signal power P^ 


(34) 


The average received signal power is given by the following equation (ref. 13): 


_ c 2 &tGrPt l_ 
lgj^Kf 2 R 2 


(35) 


where Gt an ^ Gr are th e transmitter and receiver antenna gains , P-^ is the 
transmitted power, and K is a constant which takes into account any receiver 
gain between the antenna and the input to the PLL. The table below gives the 
power gain of several representative antenna configurations. (A is the 
effective area of the antenna.) 


TABLE I.- ANTENNA POWER GAINS 


Antenna 

Power gain 

Isotropic 

Infinitesimal dipole 
Half-wave dipole 
Optimum horn 
Parabola or lens 
Broadside array 

. 1 
1.5 

1.64 

10Af t 2 / c 2 

( 6.3 to 7.5) Af t 2 /c 2 
4Af t 2 /c 2 
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Substituting equation (35) into (3*0 yields 


$ 


n 



(36) 


where Ky is equal to the quantity in parentheses in equation ( 36 ). Finally , 
substituting equation ( 36 ) into ( 33 ) results in an rms range-rate error 


cR s/Ky^ti fl + 4 f ; 2 
°R " 8 jtf t T \ 5 


^ (l - 4£ 2 )sin to d T + 1 ■ cos <*> d T 


Tll/2 

lr 


(37) 


Uncertainty in Do-ppler averaging time . - The Doppler frequency shift is 
measured by counting the number of half-cycles in a given counting interval 
T. An error in the counting time will give an error in the measured Doppler 
frequency and hence in the range rate . 


The timing pulses which determine the counting time are generated elec- 
tronically by division of the master oscillator frequency. The error in T 
is due to jitter in the frequency dividers and drift in the master oscillator. 

The frequency dividers introduce a small time delay between the applica- 
tion of the trigger pulse and the appearance of the output pulse. If the 
delay were constant, it could be compensated for; however, the effect is a 
function of transistor and circuit parameter variations, power supply voltage 
fluctuations, and extraneous signal couplings. All these factors cause a 
phase jitter to be superimposed on the timing pulse. Normally, the phase jit- 
ter is compounded in each stage of a multistage frequency divider. It is pos- 
sible, however, to design a multistage divider that has the time delay of only 
one stage (ref. 14) . Therefore, by proper design this error can be made 
negligibly small (on the order of nanoseconds). 


A more important jitter source is the master oscillator frequency 
instability. This error is analyzed by solving equation (20) for f td and 
substituting the result into equation (l) to give 


. _ cN cfo 
4f t T " 2f t 


(38) 


Equation ( 38 ) is differentiated with respect to T and the derivative 
substituted into the total differential: 



(39) 


The quantity N/2T can be obtained from equation (20) and substituted into 
( 39 ) to give 
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The first term in the parentheses of equation (40) is the range rate expressed 
in equation (l), and AT/T is the oscillator fractional frequency error whose 
rms value is the oscillator instability S s . Therefore, the rms range-rate 
error due to uncertainty in the duration of the counting interval is 

- (* + (Ul) 


It should be noted that this error is not correlated with the error of equa- 
tion (l6). The frequency variation Af/f which gives rise to equation (l6) 
is the random frequency change which occurs during the time delay t between 
the transmission and reception of the signal, whereas the Af/f which gives 
rise to equation ( 4l) is the difference between the actual frequency and the 
assumed (nominal) oscillator frequency during the count period T. These two 
errors are not the same, hence their effects on the range-rate error are 
uncorrelated. However, the rms values of both of the error contributions 
are proportional to the same fractional frequency variance, S s . 

Uncertainty in the speed of light . - It can be seen from equation (l) that 
the conversion of Doppler shift to range rate requires a knowledge of the 
velocity of light c. The velocity of light is, however, not known exactly. 
Reference 15 gives the three-sigma fractional uncertainty Ac/c as 10 “ 6 . 
Therefore, the standard deviation of the fractional uncertainty Ac/c is 
3. 33X10 - 7 . 


The error in range rate due to the uncertainty in the knowledge of the 
velocity of light can be found by differentiating equation (l) with respect to 
c and substituting the result into the total differential: 


A * dR A f td 


Ac 


c ^td Ac A Ac 


( 42 ) 


Substituting the numerical value of Ac/c into ( 42) gives the standard 
deviation in range rate due to the uncertainty in the speed of light. 

<T£ = (3-33 x 10' 7 )R (43) 

This is a bias -type error. 

Deterministic Errors 

Averaging time error *- The range rate measured by the Doppler cycle - 
counting process is an average range rate over a finite averaging time inter- 
val T. This average value is equal to the true range rate only at some 
unknown instant of time in the interval T. Usually, the average range rate 
is considered to be the true range rate at the middle of the averaging time 
interval T. Hence, there is an error 



AR = R - Ra 


(44) 


where R a is the average (measured) range rate and R is the true range rate 
at the middle of the interval T. Clearly the averaging time error is depen- 
dent on the satellite orbit geometry. An expression is derived in this sec- 
tion assuming a circular orhit. This geometry was selected because it gives 
the most insight into the problem. The more practical cases of the elliptic 
and hyperbolic orbits are derived in appendix B. 

Consider the coordinate system of 
figure 5» The coordinate axes are 
alined with the xy plane in the orbit 
plane of the satellite and with the 
beacon located at a point B above the 
y axis (x = 0) on the lunar surface. 
The moon is assumed not to rotate with 
respect to the coordinate system. The 
beacon is located by the vector 

= r B cos 3j + r B sin 3k (45) 


The satellite is located by the. vector 
r = (r Sin 9t)i + (r cos 0t)j 

(46) 

Note that when t = 0 the satellite is 
at the point of nearest approach to the 
r is found using the definition of the 

>s 3 cos 0t (47) 


Application of the law of cosines to triangle OBS and of equation (47) 
results in . 

R 2 = r B 2 + r 2 - 2r B r cos 3 cos 0t ( 48) 

The law of cosines applied to triangle OPS yields 

p 2 = r 2 + (r-g cos 3) 2 ~ 2rr B cos 3 cos 0t (49) 

Equation ( 49) is differentiated with respect to time to give 

pp = err-g cos 3 sin 9t (50) 

A relationship between p . a nd R can be obtained by use of the Pythagorean 
theorem on right triangle PBS 

R 2 = (r B sin 3) 2 + P 2 (51) 



Figure 5.- Circular orbit .geometry. ' 


beacon. The angle if between rg and 
dot product, with the result that 

cos \|r = c 
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Equation (5l) is differentiated with respect to time and solved for R. 


R = 


££ 

R 


( 52 ) 


Then equation (50) is substituted into (52) to give R as a function of range 
and time. 


0rr-r> cos (3 sin 0t 
» = ' R 


(53) 


Range. can he eliminated from (53) with the aid of equation (48). 

0rrg cos 3 sin 0t 

R — loo • 

vr„ + r - 2rr_ cos 3 cos 0t 


(54) 


B 


B 


Equation (54) expresses the instantaneous range rate of the satellite with 
respect to the beacon as a function of time. The average range rate over £ 
time interval T is found by integrating equation (54). 


T 

t+| 


T 

t+ i 


2 . 

% = T J- | Rdt= T / t _T^TF7 


Orr-r. cos P sin 0t 

B ?= dt (55) 


r 2 - 2rr T) cos 3 cos Qt 
£> 


With the change of variable u = cos 0t , equation (54) can be integrated. 5 
The result is 


Ra = 


J2rr s cos p 


T 


{[ 


r ' r B 2 + r2 ' •/ T 

2rr cos p COS + 2 
J3 x 


1/2 


[: 


r B ^ + t2 
2rr^ cos 3 


- cos 0i 


<• - m 


( 56 : 


R a = 


\ / 2rr-g cos 3 


a/ 2 


- cos 


Equation ( 56 ) can be simplified by trigonometric identities so that 

0T • QT • 

— cos 0t + sin — sin 0t ) 

\ 1 / 2 ] 

V J (57) 


p p 

t-q + r ^ 

LV 2rr B 00 s ^ 


r B 2 + r 2 


0T • ■ 

2rr>g cos P " C0S T COS 0t " sln T sin 0t 


5 Integral 191 * 01 , page 46, reference l 6 . 
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Time can be eliminated from equation (57) by solving equation (48) for 
cos 9 t. 


cos 9 t = 


2rr B cos P 


( 58 ) 


and equation (53) for sin 0t, 


RR 


sin 9t = 7 

0rr B cos P 

and substituting (58) and (59) into (57): 


(59) 


n/ 2rr^ cos |3 

y t b 2 + r2 

C\J 

1 

CM 

U 

+ 

CM 

0T 

RR . 0T\ 

T 

\2rr B cos P 

2rr B cos P 

cos 2 

5rr B cos P Sin 2 ) 


l/2 


r B 2 + 


[ B 


+ r 2 - R 2 


RR 


2rr B cos P 2rrn cos P 


B 


• 0T\ 

2 grr B cos p’ sin 2 J 


9 T 

cos — - £ 


1/2' 


(60) 


For a typical lunar orbit the angular velocity .9 is of the order of 
10 -3 rad/sec. . Therefore, 0T/2 « # 1 and equation (60) can be simplified by 
letting sin 0T/2 = 0T/2 and cos 0T/2 = 1. The result is 


R a = | (R 2 + TRR) 1/2 - I (R 2 - TRR) 1/2 


(61) 


The averaging time error for a circular orbit, using equation (44), is equal 
to 


AR = R - | (R 2 + TRR) 1/2 + I (R 2 - TRR) 1/2 


(62) 


For the special case where R < r/t (this is generally true for most 
cases of practical interest) the quantities (R 2 ± TRR) 1/2 in equation (62) 
can each be expanded by the binomial theorem. 


(R 2 ± TRR) 


dU/ 2 


- R + 1 B m _ 1 R 2 T 2 + 1 R 3 T 3 

■ R ± 2 RT 8 R ± 16 T' 


5 R 4 T 4 + _7_ r 5 t 5 
128 R 3 “ 256 R 4 


• • • 


(63) 


Substitution of (63) into (62) yields 
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(64) 


AR = - 


T 2 R 3 TT 4 ^ 5 

"SF ' 128R 4 - * * * 


If R is sufficiently less than R/T so that the second term of (64) is 
negligible, the range-rate error due to the averaging time error in a circular 
orbit becomes 



AR = - 


T 2 R 3 

8r 2 


(65) 


Figure 6.- Range rate averaging time error for 
satellite in circular lunar orbit; 


. altitude = 200 km, 0=0 

propagation time , and is detected to yield R(t). 
rate at the instant of reception is R(t + t/ 2). 
in the range -rate measurement equal to 


In figure 6 the averaging time 
error is plotted against the angular 
displacement 0 of the satellite in 
its orbit (0=0 corresponds to the 
point of closest approach to the bea- 
con) for various values of Doppler 
averaging time, T. 

Finite propagation time error . - 
The Doppler -shifted frequency trans- 
ponded by the beacon to the satellite 
is related to a range rate R(t). This 
signal arrives at the satellite at time 
t + t/2, where t/2 is the one-way 
However, the true range 
Therefore, there is an error 


AR 

= R^t + 

D 

- R(t) 




t+l 

AR 

= R(t) + 



T 

t+k 


' t 

AR 

-l 

*• 

R 

dt 


( 66 ) 

(67) 

(68) 


In the case of a lunar satellite typical propagation times will be very 
small, usually no more than 10" 2 second. Therefore, equation (68) can be 
simplified to 


“ - T = 8 (D 


(69) 


For a circular orbit equation (53) can be rearranged to give 


RR = 0rr_ cos 0 sin 0t 
B 


(70) 


and differentiated with respect to time to give 
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RR = 0 2 rr-g cos P cos 0T - R 2 


(71) 


Substituting equation (7l) into ( 69 ) and eliminating cos 0t with 
equation ( 58 ) gives the finite propagation time error for a circular orbit. 


“ - £ <■** 



. 8, deg 

Figure 7*- Range rate finite propagation time 
error for satellite in circular lunar 
orbit; altitude = 200 km, (3 = 0°. 


Figure 7 shows the finite propaga- 
tion time error in range rate versus 
the angular displacement of the satel- 
lite for a circular orbit. It is seen 
that this error is a minimum when R 
(and hence the propagation time) is a 
maximum and increases to a maximum at 
0=0 when R is a minimum. This is 
because as 0 ..decreases , the range 
acceleration R increases faster than 
the range.. R decreases, so that the 
product RR in equation ( 69 ) increases 
as 0 decreases. In this sense the 
R term of equation ( 69 ) dominates this 
error. 


RANGE 


The Ranging Principle 

If the interrogator carrier signal is phase modulated with a low- 
frequency sinusoidal tone, it is possible to measure range. This is possible 
because the signal transmitted by the interrogator requires a finite t im e t 
to complete the trip to the beacon and back to the interrogator. Therefore, 
there will be a phase difference cp m between the transmitted signal and the 
received signal equal to ~ 

Tm ~ ^m 1 ” (73) 


The two-way propagation time t 
equation (73) becomes 


is approximated by t = 2R/c 


c Pm = 


4jtf m R 


so that 


(74) 


Equation (7*0 can be solved for the range to give the range equation 


R . 


c <P m 

4jtfm 


(75) 
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The maximum unambiguous phase angle cp m is 2jt radians, which results in 
a maximum unambiguous range of 


r _ _ c _ ^m 

^max kjtf m 2f m 2 


( 76 ) 


Therefore, the wavelength of the ranging tone should be equal to or greater 
than twice the maximum expected range if there is no way of resolving the 
range ambiguity. For example, an 800 Hz ranging tone has a maximum 
unambiguous range of 187.3 km. 

However, the phase angle cp m cannot be measured exactly. For the par- 
ticular system described in reference 1 the overall phase measurement accuracy 
is about one percent of full range or 3*6° in phase. This means that the 
range resolution of the 800 Hz tone would be 1.873 km. To improve the resolu- 
tion, the carrier is phase modulated with additional higher frequency ranging 
tones. The low-frequency tones are used only to remove the range ambiguity of 
the high-frequency tones. For example, a 100 kHz high-frequency ranging tone 
would have a resolution of 15 m. 


Description of the Ranging System 

The basic block diagram of the interrogator portion of the ranging system 
is shown in figure 8. The portion within the dashed line is the same 
equipment as that enclosed within the dashed line of the range rate system in 
figure 1. 



Figure 8. - Ranging system. 
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The ranging tones are generated hy frequency division of a stable master 
oscillator frequency. The output of the frequency dividers is several coher- 
ent. square wave ranging tones (two, f mi and f m? , are shown in fig. 8) which 
are filtered to provide sine waves and fed to individual phase shifters which 
adjust the phase of each tone and remove any differential phase shift that may 
be present in the system. The properly phased coherent ranging tones then 
phase modulate the carrier and are transmitted to the transponder. The trans- 
ponder receives the signal, amplifies it, and retransmits it back to the 
interrogator. The beacon transponder equipment is the same as that used in 
the range-rate measuring system. 

The receiver PLL within the dashed portion serves a dual role as a nar- 
rowband tracking filter and as the first IF stage. Therefore, the PLL 
output is a carrier at a constant IF frequency and phase modulated with the 
ranging tones. These are passed through a bandpass filter bank to separate 
the individual ranging tones , and each tone is then tracked with a second PLL 
to further improve the signal-to -noise ratio. It should be noted that each 
ranging tone has a slight Doppler shift. This Doppler shift could be removed 
by a phase -locked loop, but it is a deterministic error and so can more easily 
be removed by a computer during data reduction. This error is analyzed in a 
later section. 

The phase angles of the ranging tones are measured in either of two ways: 

(1) by comparison with a reference (transmitted) tone in a phase detector, or 

(2) by means of an all digital pulse counting technique. In the first case, 
the phase detector output is an analog voltage proportional to the phase dif- 
ference between the transmitted and received signals. In the second case, the 
received ranging tones, by means of a stop-pulse generator, determine the 
length of time that a counter gate is open. This time interval determines the 
number of pulses of frequency f c which are counted by the pulse counter, and 
this count is proportional to the range. A detailed description of the 
operation of the stop-pulse generator is given in references 1 and 2. 


Stochastic Errors 

Master oscillator short-term instability . - The ranging system works on 
the principle of measuring the phase difference between the transmitted and 
received ranging tones. Clearly then, the ranging frequencies must have good 
phase, and hence frequency, stability, if an accurate phase comparison is to 
be made. 

Consider the transmitted and received ranging tones at the input of the 
phase detector. The reference (transmitted) ranging tone can be represented 
as 


yt(t) = sintwjjjt + cp(t)] (77) 

where is the ranging tone radian frequency and cp( -fc) is the phase noise 

of the oscillator output. The received tone is 


2k 



y r (t) = sintw^t - t) + q>(t - t)] 


(78) 


The output of the phase detector is proportional to the phase difference 


cPm(t) = WmT + cp(t) - <p(t - t) 


(79) 


The range, from equation ( 75) > is 


R = + 2 ^ Mt) " c P( t " (80) 

The rms error in the range measurement due to oscillator instability is 
therefore equal to 

cr R = ^ jjpltp - ap(t)<p(t - t) + cp(t - t) 2 J (8l) 

where the bar means "average." The technique used in the discussion of range 
rate short-term oscillator instability can be used to write equation (8l) in 
terms of the autocorrelation function of the phase noise cp(t) as 


d R = -pP- [R( 0 ) - R(t )] 1/2 (82) 

R J&m 


Over the narrow bandwidth of the ranging tone receiver the phase fluctua- 
tion cp(t) is approximately white noise. Therefore, assume it has a two-sided 
power spectral density of low pass filtered white noise 


V u) 


U ) 2 + V 2 


(83) 


The autocorrelation function ft(r) can be found by the inverse Fourier 
transform of Gp(w) 


ft(-r) 



(84) 


Because of even symmetry 


R(t) 


?£i 

at 



COS T 
CO 2 + v 2 


(85) 
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Equation (85) can be integrated with tables 6 to yield 


R(t) 



(86) 


A physical interpretation of v can be found by considering the noise 
equivalent bandwidth (ref. 17) , B: 


Total power 
Peak power 



(87) 


The noise equivalent bandwidth is the bandwidth of a rectangular power spec- 
tral density with height G^O) and with an area equal to that under Qp(oj) . 
Substituting equation (83) into (87) yields 


r 2 ° Slt „ ■ rad - v tt 
B = 0 o /,, = Vrt = rHz 

2cr^/v sec 2 


( 88 ) 


Thus, v is equal. to twice the noise equivalent bandwidth of the phase 
fluctuation. 

The physical interpretation of Op 2 is easily seen by noting that 
ft(0) = cTp 2 = mean square power. Since the phase fluctuation has zero mean, a p 2 
is the variance of the phase fluctuation. 

Equation (86) is substituted into (82) to obtain the rms range error 
due to short-term oscillator instability. 


R 




(89) 


Most oscillator phase noise is approximately white noise; therefore, v is 
large. In this case equation (89) can be simplified to 


cg p 


(90) 


If the oscillator rms fractional frequency instability is S g = Af^j/fjjj, then 
the rms phase deviation ffp is 

<»> 

6 Integral 859. 001, page 224, reference l6. 
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Substituting equation (91) into (90) and letting t = 2R/c yields 

<X r =n/ 2S s R (92) 

This is a noise -type error. 


Master oscillator long-term instability .- The error in range due to the 
master oscillator long-term instability is a bias -type error caused by an 
uncertainty in the modulating frequency f m . The error is calculated by 
differentiating equation (75) with respect to f m and substituting the 
derivative into the total differential, 


AR = ^Af m 
UI m 



(93) 


Since the first term in parentheses is the range R and the second term is 
the long-term rms fractional frequency stability factor equation (93) 
becomes 


°E - s l R ^b) 

Quantization error .- In the. digital range measuring scheme the master 
oscillator generates a start pulse which opens an electronic gate and allows a 
pulse counter to count pulses of frequency f c . At a later time, equal to the 
two-way propagation time t, the coherent ranging tones generate a stop pulse 
which closes the gate and stops the counter. The counter counts zero cross- 
ings of the pulses which pass through the gate; therefore, the counter has a 
resolution of only half a cycle. Since the stop pulse will generally not 
occur exactly at a zero crossing, the counter reading will always be a frac- 
tion of a half -cycle less than the true count (i.e., the counter counts to the 
nearest integer less than the true count) . 

In the time r the counter (initially cleared) will count to 


N - rf c = ^ f c (95) 

The error in range due to a counting error AN is found by solving equa- 
tion (95) for range, differentiating with respect to N and substituting the 
derivative into the total differential: 

. 4 R= i 4 ,= i a < 96 > 
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There is a maximum count error AN may of one count (negative) at the end of 
the count period. Therefore, the maximum range error, from equation (96), is 
ARmax = _ c/2f c . There is no count error at the beginning of the count period 
because the start pulse is coherent with the transmitted ranging tones. 

The method of analysis is analogous to the range-rate quantization error 
analysis. Thus, the mean range quantization error can be shown to be AR ^y /2 
and the standard deviation to be AR mx /2 -J3. Therefore, the mean range 
quantization error is 

mean (AR) - (97) 


and the standard deviation of the range quantization error is 


^ = 4 71 f c 


(98) 


To obtain zero mean errors in range measurement, the mean error (eq. (97)) 
should be added to the range measurement. This could be done in the computer 
by adding half a count to N in equation (95) before computing the measured 
range . 

Phase -locked loop error . - Since the ranging principle depends directly on. 
the measurement of a phase angle difference, imperfect phase tracking of the 
ranging tone signal by the ranging tone PLL will cause an error in the range 
measurement. This section considers this error. The additive noise at the 
input of the PLL is assumed to be Gaussian and white. The PLL is assumed 
to be a second-order active loop similar to that discussed in the range-rate 
section. 

The reference (transmitted) signal tone can be represented as 

E t (t) = sin w m t " (99) 


and the received signal at the same instant is 

E r (t) = sin w^(t - t) + N(t) (lOO) 


where N(t) represents the additive Gaussian white noise. Jaffe and Rechtin 
(ref. 18) have shown that when the PLL is linearized; the additive noise can 
he represented as a phase noise cp(t) if the phase error is small. Therefore; 
the received signal can he represented as 
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E r (t) = sin[c%(t - t) + cp(t)] 


( 101 ) 


The reference and received signals are phase detected to yield a signal 
proportional to the phase difference 


<P m (t) = WmT - cp(t) 


( 102 ) 


Substituting equation (102) into equation (75) yields 


ct _ ccp(t) 

2 2^m 

The second term of equation (103) represents the error in the range 
measurement. The rms range error is 


(103) 






(104) 


where the bar means "average," and ft( 0) is the autocorrelation function ft(t) 
with the correlation time equal to zero. But fl(0) in equation ( 104) is 
identical in form to ft(0) in equation (32). Therefore, the rms range error 
due to additive phase noise is 


c /*R M n l+J^fY /2 

a R Sn\ 2 5 J 


(105) 


If Nr is the power spectral density of the phase noise cp(t), then letting 
<&R = KrR 2 , as i n the case of the range-rate PLL error equation (36), 
equation ( 105 ) becomes 


cR / ¥n 1 + 4C 2 Y /2 

" 8icf m V 2 5 J 


(10 6 ) 


Phase detector error .- A phase detector does not measure a phase differ- 
ence without error. Therefore , if the range is determined by measuring the 
phase difference between the received and transmitted ranging tones > the phase 
detector error must be considered. 

The error in range produced by the inaccuracy of the phase detector is 
found by differentiating equation (75) with respect to and substituting 

the derivative into the total differential. 

^ = l§5 ^ < 10 i> 
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where Acp m is the amount of phase error introduced by the phase detector 
expressed in degrees. If 0 p D is the rms value of Ap m , the rms range 
error due to the phase detector error is 


= cg H) 
a R 720 f m 


( 108 ) 


This error exists in a ranging system only if range is measured by the phase 
detector method. If range is measured with the digital pulse -counting method, 
the system will have quantization error instead. 

Calibration drift error . - The ranging system requires the correct phase 
relationship between the ranging tones. Adjustable phase shifters are, there- 
fore, incorporated in the circuitry to insure the proper phase relationships 
and to remove the phase bias errors inherent in the transmitting and receiving 
equipment. A range measurement error results whenever the phase shift 
calibration drifts. 

Calibration drift is a significant source of error with ranging equipment 
subjected to environmental extremes. Such errors vary widely with equipment 
and environment, and usually can be determined only experimentally. A repre- 
sentative 3 d value for the overall system calibration error (ref. 3 ) is 
0.869°. Thus, the standard deviation in phase drift cJqq is O. 289 0 . The 
standard deviation in range is found by substituting the rms phase drift 
0 QD into equation (107). The result is 


= c(J CD 
a R 720f m 


(109) 


This is a bias-type error (provided that the phase drift does not change 
significantly during the time period of interest). 

Uncertainty in the speed of light .- It can be seen from equation (75) 
that the conversion of measured phase angle cp m to range requires a knowledge 
of the velocity of light. However, since there is an uncertainty in the exact 
value of c, there is a range error equal to 


tew 


cc Pm Ac _ Ac 
c " R T 


(no) 


where 8 r/ 8 c is found by differentiating equation (75). The la uncertainty 
is 3 . 33 x 10 -7 , so that the rms range error due to the uncertainty in the 
speed of light is 


d R = (3.33><10- 7 )R 


This is a bias -type error. 
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Deterministic Errors 


Finite •propagation time error . - There exists a finite propagation time 
error in the range measurement analogous to that described for the range rate. 
(For the following analysis refer to figure 5*) 

Equation (48) gives the beacon -to-satellite range at the instant the 
beacon transponds. The satellite measures this value which will be called R^ 


R 


m 


-J 


r 2 + 
r B + 


r 2 - 2rrg cos P cos 9t 


( 112 ) 


A better estimate of the range is 


R = Jrg 2 + r 2 - 2rr B cos 3 cos §(t + R m /c) 


(113) 


where t + R m /c is the time at which the ranging tones are received. 
Normally 6>R m /c <<: 1 so that equation (113) simplifies to 


A 

R 



- 2rrg ( cos 9 1 



sin 0t 


cos p 


(114) 


Equation (112) can be substituted into (ll4) to give 


R 



2rrg0 cos P 
: — R m sin 9t 


(115) 


The finite propagation time error is equal to 


AR = R - R m 


(116) 


Therefore , 


AR = 



2rrg0 cos P 


R m sin 0t - R, 


m 


(117) 


The explicit dependence of AR on time. is removed by substituting 
equation (59) into ( 117 ) with the .result that 



(118) 
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Since 2R/c « 1 for typical orbits, the square root term of equation (ll8) can 
be expanded with the binomial expansion 



2R • 1 

c 


R 

c 


(119) 


Therefore , 

AR = 55 (120) 



Figure 9.- Range finite propagation time error 
for satellite in circular lunar orbit; 

3 - 0 °. 


Figure 9 shows the finite propagation 
time error in range versus the angular 
displacement 9 of the satellite for 
various circular orbit altitutdes. The 
curves -terminate at values of 9 at 
which the satellite recedes below the 
lunar horizon. 

Doppler shift of the ranging 
tones . - If r is the two-way propaga- 
tion time, the phase angle of a ranging 
tone with no Doppler shift is 

cp m = 2itf m T (121) 

For the Doppler shifted signal the 
phase angle is 


cPmH - ^(^m + ^md) T (122) 

The error in the phase angle measurement is 

Acpm = Tmd ” Tm = ^^md T (123) 


Substituting equation (123) into (75) and letting r = 2R/c gives 



But from appendix A 


^m + ^md 


c_+R 
c - ft 


L m 


(124) 


(125) 
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which can he solved for fmd/ fm to give 


f md _ 2R 
f m c - ft 


(126) 


Therefore, the error in range due to a Doppler shift in the ranging tone is 
equal to 


AR = 



TOTAL ERROR MODEL NUMERICAL EXAMPLES 


(127) 


This section gives expressions for the root -sum-square (rss) range and 
range-rate stochastic errors for typical system parameters. The deterministic 
errors are not included because they can be removed during data reduction. If 
the individual stochastic errors described in the previous sections are con- 
sidered to be (lb) standard deviations of independent random processes, then 
the overall range and range -rate variances can be found by summing the 
individual variances . 


The following typical parameters are used for the numerical example: 


fc 


G r 

<*t 

K 

W R 

*t 

Sl 

S s 

T 

t 

a CD 

a PD 


10MHz 

100kHz 

200kHz 

1.7GHz 

1 

1 

1 

2.5X10" 21 w/Hz 

4X10 -2 ° W/Hz 

1 W 

1X10 " 6 

1X10 " 9 

1 sec 

0.5 

0 . 289 ° 

3° 

6.28 rad/sec 
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Assume that the satellite is in a circular orbit about the moon with an 
altitude of 200 km, and the beacon latitude is 0° (overhead pass). 


The rss range-rate error is: 


r (3X10 _lo )R for 0 < R < I.5XIO 8 

R t 4. 5X10 " 2 for R > 1.5X10 8 

Master oscillator short - 
. term frequency instability 
(eq. (16)) 

+ (1X10 _12 )R 2 

Master oscillator long-term 
frequency instability 
(eq. (19)) 

+ 3. 24x10 " 4 

Quantization (eq. (25)) 

+ (24.4xi0 _2O )R 2 

Phase -locked loop (eq. (37)) 

+ (R + 1.77X10 4 ) 2 10" 18 
31/2 

+ (1.11X10 _13 )R 2 ^ m/sec 

Uncertainty in Doppler 
averaging time (eq. (4l)) 

Uncertainty in speed of 
light (eq. (43)) 


(128) 

where R is in meters and R is in m/sec. It 
in range rate is caused by the master oscillator 
instability and by the quantization process. 

is seen that the largest error 
short-term frequency 

The rss range error is 


cTr = -|(2xlO" ia )R 2 

Master oscillator short- 
term frequency instability 
(eq. (92)) 

+ (1X10' 12 )R 2 

Master oscillator long-term 


frequency instability 
(eq. (94)) 


+ (2.28xi0~ 12 )R 2 


l8.7 (quantization) 

17.4 (phase detector) 

+ 1.45 

] 1/ 2 

+ ( 1. 11x10 ~ 13 ) R 2 f meters 


Phase -locked loop 
(eq. (106)) 

(eq. (98)) 

(eq. (108)) 

Calibration drift 
(eq. (109)) 

Uncertainty in speed of 
light (eq. (ill)) 

(129) 
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It is seen that the range error due to oscillator short-term frequency 
instability is negligible compared to the other error sources. 

Equations ( 128 ) and (129) are plotted against angular displacement 6 in 
figures 10 and 11. For the example orbits the range R is always less than 
1.5X10® m, so that the first term of (128) is dependent on R. 



Figure 10 . - Total range rate error for satellite 
in circular lunar orbit; 3 = 0°, T = 1 sec. 



Figure 11.- Total range error for satellite in 
circular lunar orbit; 3=0°. 


CONCLUDING REMARKS 


The descriptions of tracking system errors given constitute a reasonably 
complete error model for an interrogator-transponder type range and range -rate 
system. No attempt has been made here to evaluate the relative significance 
. of the various types of errors since this depends upon the particular system 
and tracking mission considered. Generally, when the results are applied to a 
specific case, it will be found that only a few of the errors need to be 
included. 


Ames Research Center 

National Aeronautics and Space Administration 

Moffett Field, Calif., 9^035, Sept. 6 , 1967 
125 - 17 - 05 - 01 - 00-21 
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APPENDIX A 


DERIVATION OF THE DOPPLER SHIFT EQUATION 


The satellite approaches the beacon transponder and transmits with 
frequency f t - The transponder receives this signal at frequency 


= (c + R 
\c - b 


vl/2 


(Al) 


in accordance with the standard equation describing the Doppler effect 
(ref. 19). The beacon retransmits with frequency 


W 


f t + f L 


(A2) 


Finally, the satellite receives the transponded signal at frequency 




+ fr 


c + R n1/2 


'J \c 


(A3) 




L 2 VC - 


c - 


(A4) 


By definition 


+ 2f td “ lVTTJ f t 


(A5) 


and 


fL + f„ 


c + r\ 1/2 . 


L Ld \c - iy L 
Equation (A5) is solved for the range rate R giving 

• cf td 


R = 


f td + f t 


(A6) 


(A7) 


However, since f^ » f^ equation (A7) can be simplified to 


36 



(A8) 


R « 


cf td 


If the vehicle and beacon are moving away from each other, the result is 
exactly the same. 
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APPENDIX B 


AVERAGING TIME AND FINITE PROPAGATION TIME ERRORS 
FOR ELLIPTIC AND HYPERBOLIC ORBITS 


Figure 12 illustrates the orbit configuration to be considered. The 
rectangular coordinate system has its origin located at one focus of the conic 



Figure 12. - Orbit configuration. 


and is arranged with its xy plane in the satellite's orbital plane and with 
the y axis in the direction of the perilune. The beacon is located by 
longitude A and latitude (3. The satellite is located in the orbit by the 
true anomaly 0. The elliptic orbit will be considered first. 

Three basic equations from the celestial mechanics of elliptic orbits 
that are useful in the subsequent analysis are given below without proof . 1 


1 These are, respectively, equations 3 - 35 , 3 - 36 , and 3-40 of reference 20. 
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r cos Q = a(cos E - 

e) 

(Bl) 

r sin 9 = a vl - e 2 

sin E 

(B2) 

nt = E - e sin 

E . 

(B3) 


where a is the semimajor axis of the satellite orbit , e is the e ccentr icity, 
E is the eccentric anomaly, and n is the mean angular velocity, v n/a 3 . 

Referring to figure 12 it is seen that: 


r- = r,, cos P sin Ai + r^. cos P cos Aj + r^, sin Pk 
n d si si 


r = r sin 0i + r cos 9 j 


and 


R = r - r 
B 


(B4) 

(B5) 

(b6) 


Substituting equations (B4) and (B5) into (b 6) yields 


R = (r sin 9 - r^ cos P sin A)i + (r cos 9 - r B cos P cos A) j - r^ sin Pk 

(B7) 

Equation (B7) is written in terms of the eccentric anomaly instead of the true 
anomaly by substituting equations (Bl) and (B2) into (B7) • The result is 

R = ^a \/l - e 2 sin E - r B cos P sin A^i 

+ J"a(cos E - e) - r B cos P cos A^j - r B sin Pk 


(B8) 


Equation (B8) is differentiated with respect to time to yield the range rate 


R = a '/l _ 


dE t 


dE - 


e cos E — i - a sin E — j 
AA - dt 


dt 


(B9) 


ass uming that the Moon does not rotate with respect to the coordinate system. 
The derivative dE/dt is obtained from equation (B3) and substituted into 
equation (B9) so that the range rate becomes 
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a -an sm E t . an n/i - e 2 cos E — 

R = ^ x + = — j 

1 - e cos E 1 - e cos E 


(BIO) 


The average range rate over an averaging time T is found by evaluating the 
integral 


T 

— 1 r ^*2 V T P . “ P 

* a = T J T B dt = T J T * dt + T J T * dt ^ B11 ^ 

■c 2 t_ 2 t_ 2 

Substituting the components of (B9) into (Bll) and integrating yields 

(B12) 


T 

* 


T 
‘ t+7T 


R a - 


n/1 - 
T 


(sin E 2 - sin Ei)i + ^ (cos E 2 - cos Ei)j 


where Ei is the eccentric anomaly corresponding to t - (t/ 2); E 2 is the 
eccentric anomaly corresponding to t + (t/ 2); and Ei and E 2 are found by 
solving equation (B3) • However, one cannot solve equation (B3) explicitly for 
E; therefore, some numerical scheme such as the Newt on -Raphs on method must be 
used. 


An approximate method giving the average range rate R a as an explicit 
function of the true anomaly Q can be obtained by differentiating 
equation (B3) and substituting the derivative into the total differential 


AT1 dE n At 

AE — ^ At i — 

dt 1 - e cos E 


(B13) 


Then letting 


At = T/2, the approximate expressions for Ei and E 2 are 


Ei 


= E nT/2 

1 - e cos E 


E2 


E + 


nT/2 

1 - e cos E 


(B14) 

(B15) 


Equations (Bl4) and (B15) are substituted into equation (B12) and the result 
is simplified "by trigonometric identities to yield 


T 2a -JT 


R a = 


R 


si / nT / 2 A cs El - f 81 /_SSA_A 

\1 - e cos E \J T \1 - e cos E \J 


(sin Ej 
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(Bl6) 



Finally, equation (Bl6) is expressed in terms of the true anomaly Q by using 
the following equation 2 relating 6 with the eccentric anomaly E: 


E = 2 arc tan 



(B17) 


Now the averaging time error can be calculated using the equation 


AR = R - 


(B18) 


In terms of the vector components of R and R a , equation (Bl8) can be written 
as 

AR = n/x 2 + Y 2 - -/X a 2 + Y a 2 (B19) 

• • • • 

where X, Y, Xg., Y a are available from equations (BIO) and (Bl6) . 


The range and range-rate finite propagation time errors are given approx- 
imately by equations (120) and (66), respectively. Evaluation of these equa- 
tions requires that the magnitudes R, R, and R be known. The range 
magnitude R can be obtained from equation (b 8) , and the range -rate magnitude 
R can be obtained from equation (BIO) . An expression for the range accelera- 
tion R can be obtained by differentiating equation (BIO) with respect to 
time and eliminating dE/dt with the aid of equation (B3) . The result is 



an 2 \/l - e 2 sin E — _ an 2 (e - cos E) — 
(l - e cos E) 3 1 (l - e cos E) 3 J 


(B20) 


The analysis of a hyperbolic orbit is analogous to that of the elliptic 
orbit; therefore, only the results are stated here. The derivation of the 
results makes use of the following three equations from the celestial 
mechanics of hyperbolic orbits (note the similarity to eqs. (Bl), (B2) , and 
(B3)): 

r cos 9 = a(e - cosh F) (B2l) 

r sin 6 = a -J e 2 - 1 sinh F (B22) 


nt = e sinh F - F (B23) 


where F is the eccentric anomaly in a hyperbolic orbit. The results are 


2 For derivation see page 51 of reference 20. 


kl 


R = (a v e 2 - 1 sinh F - r-g cos P sin 7\) i 


+ [a(e - 

cosh F) - 

• r B cos p 

cos 

A] j - r B sin Pk 

(B24) 

i an 4 e 2 - 1 

cosh F t . 

_ 

an sinh F — 

( "ROR 

n = 

e cosh F 

- 1 1 + 

e 

cosh F - 1 J 


-an 

v/e 2 - 1 sinh F t . 

an 2 

: ( cosh F - e) — 

( "R \ 


5 cosh F - 

.1)3 ^ 

(e 

cosh F - l) 3 J 


a **/e 2 - 1 
T 

(sinh F 2 

- sinh Fi) 

i - 

• ^ (cosh F 2 - cosh Fi)j 

(B27) 


T 

where Fi is the eccentric anomaly corresponding to t - -x; F 2 is the eccen 

T ^ 

trie anomaly corresponding to t + ^ an <i F i and F 2 are found hy solving 
equation (B23) or by using the following approximate relationships: 


Fi = 


nT/2 

e cosh F - 1 


(B28) 


F 2 = F + 


nT/2 

e cosh F - 1 


(B29) 


These result in an approximate expression for the average range rate equal to: 
slnh (e collhr - l) COsh ** + f slnh (e col!hF - l) sini 

(B30) 

Finally, the true anomaly 9 is related to the eccentric anomaly F by the 
following equation: 

vl/2 


F = 2 arc tanh 


( e - lY /2 e\ 

Lvttv 2J 


(B 31 ) 


The averaging time error and range and range -rate finite propagation time 
errors are found using the above equations and equations (B19 ) > (120 ) , and 
( 66 ). 
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